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1.  Introduction 


Many  electrical  (and  for  that  matter  nonelectrical)  devices,  including  mi¬ 
crowave  sources  [1],  lasers  (2].  solid-state  mtcrowave  oscillators,  and  elec¬ 
tronic  circuits  (3],  exhibit  chaos  when  pushed  to  high  power  levels.  The 
presence  of  chaos  in  these  devices  can  limit  their  efficiency  and  usefulness. 

Recently,  several  experimental  groups  have  demonstrated  that  it  is  possible 
to  control  chaos  in  physical  devices  [2]  with  tiny  perturbations  using  con¬ 
trol  methods  that  are  based  on  the  fundamental  technique  proposed  by  Ott, 
Grebogi,  and  Yorke  [4).  A  recent  Army  Research  Laboratory  (ARL)  techni¬ 
cal  report  (5]  discusses  the  use  of  tiny-perturbation  control  for  the  transmis¬ 
sion  of  information  using  chaotic  dynamics.  The  information  is  embedded 
in  the  symbolic  dynamics  of  the  controlled  chaotic  trajectory.  (The  original 
idea  for  controlling  chaos  using  tiny  perturbations  involved  the  stabilization 
of  the  periodic  orbits  that  are  embedded  in  the  attractor  of  a  chaotic  dy¬ 
namical  system.) 

In  this  report  I  develop  a  technique  for  controlling  the  symbolic  dynamics 
[6]  of  chaos  for  the  synthesis  of  time-limited  waveforms  (pulses)  and  for 
the  synthesis  of  power  spectra  with  desired  properties.  This  technique  is  in¬ 
tended  for  the  formation  of  pulses  in  high-power  devices  using  tiny  micro¬ 
electronic  control  circuitry,  so  the  power  generation  device  is  extremely 
simple  and  efficient,  while  the  control  device  is  compact,  low-powered,  and 
fast. 

The  technique  for  forming  pulses,  like  the  technique  for  information  trans¬ 
mission.  relies  on  the  concept  of  controlling  the  symbolic  dynamics  of  a 
chaotic  system.  (One  can  view  the  control  of  periodic  orbits  as  control  of 
orbits  with  cyclic  symbolic  dynamics,  the  transmission  of  information  as 
control  of  orbits  with  aperiodic  and  stationary  symbolic  dynamics,  and 
pulse  formation  as  control  of  orbits  with  transient  symbolic  dynamics.) 

In  this  report  the  concept  of  spectral  shaping  in  a  chaotic  dynamical  system 
is  also  introduced.  The  obtainable  spectrum  is  constrained  by  the  dynamics 
of  the  system  (the  tiny-perturbation  control  does  not  alter  the  basic  topo¬ 
logical  structure  of  the  attractor),  and  is  related  to  pulse  formation.  It  is  pos¬ 
sible  to  alter  the  spectrum  of  a  chaotic  system  with  ergodic  and  stationary 
dynamics,  but  in  this  report  I  limit  the  discussion  of  spectral  shaping  to 
pulse  waveforms.  Both  ideas  introduced  here  therefore  concern  the  use  of 
chaos  to  generate  a  signal  class  that  has  not  been  addressed  previously  in 
the  context  of  chaos  control:  the  transient  time-limited  signal.  In  addition  to 
the  time  sequence  pulses,  and  the  spectra,  I  also  include  graphs  of  sinusoi¬ 
dal  waveforms  that  have  an  amplitude  modulation  described  by  a  linearly 
interpolated  continuous-time  version  of  the  pulse  waveforms.  These  graphs 
are  intended  to  illustrate  the  appearance  of  intensity  (amplitude)  modulation 
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chaos  and  controlled  chaotic  pulses  cf  intensity  modulation.  Microwave 
sources  and  lasers,  among  other  devices,  are  known  to  exhibit  this  type  of 
chaos.  (The  band-limited  chaos  that  appears  in  these  devices  is  more  likely 
to  be  ootii  phase  and  amplitude  modulation.) 

Because  different  dynamical  systems  have  attractors  with  different  topo¬ 
logical  structures  in  state  space,  the  exact  description  of  a  chaotic  system 
depends  strongly  on  the  specific  system  under  study.  Useful  descriptions  of 
chaotic  dynamical  systems  for  control  almost  without  exception  rely  on  the 
extraction  from  measured  data  of  a  simple  return  map  (a  discrete-time  de¬ 
scription  of  the  dynamics)  for  a  Poincare  surface  of  section  in  state  space.  1 
therefore  limit  the  discussion  here  to  a  specific  discrete-time  dynamical  sys¬ 
tem — the  shift  map  (7],  (The  relationship  between  a  symbolic  shift  and  a 
physical  system  has  already  been  addressed,  and  an  experiment  is  planned 
to  demonstrate  inforrmtion  transmission  using  this  relationship  [5).  This  re¬ 
port  therefore  concentrates  on  the  symbolic  description.  The  shift  map  itself 
approximately  desenbes  the  Loren?  [8]  three-dimensionai  continuous-time 
dynamical  system  for  the  commonly  used  parameter  values,  and  its  basic 
properties  appear  in  many  other  systems.) 

The  so-called  double-scroll  electrical  oscillator  that  we  have  constructed  in 
our  laboratory,  although  a  nonhyperbolic  system  of  two  collided  Rossler 
attractors,  can  be  described  by  one-dimensional  mapping  that  is  similar  to 
the  shift  in  the  ways  that  are  important  for  this  discussion.  More  complex 
systems  require  the  derivation  of  a  symbolic  system  description  from  meas¬ 
urements;  again,  this  technique  has  already  been  addressed  [5],  More  com¬ 
plex  systems,  including  die  double-scroll  system,  should  provide  even  more 
flexibility  in  waveform  synthesis. 
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2.  Background 


The  shift  map  dynamical  system  [7]  is  described  by  the  equation  x^  -  lxt 
mod  1.  It  can  therefore  be  described  as  follows  in  terms  of  a  binary  sym¬ 
bolic  dynamics.  If  the  current  state  of  the  map  is  xt  =  0 ,bfi7by..  in  a  binary 
fraction  decimal  representation,*  then  the  shift  is  alternatively  described  by 
a  shift-and-truncate  operation  on  the  binary  fraction  representation  of  the 
state  point.  To  obtain  x^.  simply  shift  the  bits  in  the  binary  fraction  for  x , 
to  the  left  bv  one  place,  and  discard  tne  bit  that  ends  up  to  the  left  of  the 
decimal  point.  Thus,  =  0 .  The  equation  =  2x„  modi  can 
thus  be  replaced  by  the  alternative  model  of  a  simple  bit  shift  operating  on 
an  infinite  binary  string.  (The  relation  between  a  symbolic  dynamical  sys¬ 
tem  and  a  physical  electrical  oscillator  dynamical  system  has  been  demon¬ 
strated  in  the  context  of  information  transmission  [5],  so  I  concentrate  on 
the  symbolic  system  here.) 

This  alternative  symbolic  description  is  appealing  both  from  an  intuitive 
and  an  analytical  standpoint:  It  is  easy  to  visualize  and  easy  to  construct  ini¬ 
tial  points  that  yield  a  desired  dynamics.  If  the  system  state  point  falls  in  the 
interval  [0,1/2),  the  system  is  sai"’  to  generate  a  binary  symbol  0,  and  if  the 
state  point  falls  in  [1/2,1],  the  system  is  said  to  generate  a  1.  Now  the  first 
bit  ;n  the  binary  fraction  for  the  state  point  determines  which  of  these  inter¬ 
vals  the  state  point  falls  in,  and  t*'e  action  of  the  shift  is  to  move  all  the  bits 
to  the  left  one  space  for  each  iteration  of  the  map.  The  symbol  sequence 
generated  by  this  map,  in  binary  fraction  notation,  is  therefore  precisely  the 
same  as  the  binary  fraction  representation  of  the  system  state  point.  Thus, 
denoting  the  system  symbolic  state  by  r ,  the  symbolic  state  is  related  to  the 
state  space  state  by  the  identity  function  One  can  therefore  specify 

the  symbolic  dynamics  by  simply  setting  the  bits  in  the  binary  fraction  for 
the  initial  point  to  the  desired  symbolic  dynamics. 

In  other  systems,  however,  the  system  state  and  the  symbolic  state  will  be 
related  by  a  function  other  than  the  identity.  (In  the  Lorenz  system,  for  ex¬ 
ample,  r(x)  is  practically  continuous  and  monotonic,  but  is  not  the  identity. 
It  can,  however,  be  transformed  into  the  identity  using  a  continuous  mono¬ 
tonic  coordinate  change.)  The  problem  of  specifying  a  time-limited  wave¬ 
form  is  now  reduced  to  the  task  of  specifying  the  bits  in  ihe  binary  fraction 
for  the  initial  point  x0.  This  conceptual  abstraction  is  justified  in  view  of  the 
fact  that  with  tiny  controls  the  symbol  sequence  that  evolves  from  the  sys¬ 
tem  can  be  viewed  as  emerging  from  the  bits  that  are  altered  below  the 
threshold  of  observability. 


'‘Briefly,  if  ihe  i'h  place  behind  the  decimal  point  is  given  by  the  binary  digit  (bit)  br  then  .he  state  point  of  the  map 
is  z  =  X“_  j  b-  2~‘ .  Each  real  number  is  thus  identified  with  an  infinite  binary  string,  and  vice  versa. 
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Fin?,  for  completeness  and  comparison,  I  briefly  describe  the  natural  dy¬ 
namics  and  spectrum  of  the  shift  map.  The  natural  invariant  probability 
density  for  the  shift  [7]  is  simply  the  uniform  density  on  the  interval.  The 
correlation  function  for  the  shift  map  signal  is  computed  using  the  formula 

C(.m)  =  j\p{x)fm(x)  dx  - 

where  fm(x)  is  the  m1*  iterate  of  the  map,  and  £[jr]5  =  1/4  is  the  squared 
mean  value  of  the  signal.  The  correlation  function  is  computed  to  be  C(rrt)  = 
(1/12)2-|J— .  (The  correlation  function  evaluated  for  m  =  0  is  equal  to  1/12, 
the  variance  of  a  uniformly  distributed  random  variable  on  the  unit  inter¬ 
val.)  The  unnormalized  power  spectral  density  is  the  Fourier  transform  of 
the  autocorrelation  function, 

S(e>)=  j^COn)*"-". 

The  unnormalized  power  spectral  density  is  given  by  S(co)  = 
(l/4)/(5  -  4  cos  co).  This  power  spectral  density,  because  it  is  for  a  discrete¬ 
time  signal  (a  sequence  of  numbers),  is  periodic  with  period  2k.  Note  that 
the  power  density  for  this  signal,  because  it  is  not  delta  correlated,  is  mere 
band  limited  than  that  of  a  uniform  random  variable,  which  has  a  flat  spec¬ 
trum  because  it  is  delta  correlated.  Figure  1(a)  shows  a  time  sequence  of 
100  points  generated  by  the  natural  dynamics  of  the  shift.  The  power  spec¬ 
tral  density,  normalized  so  that  the  integrated  spectral  power  is  unity,  is 
shown  in  figure  1(b).  A  waveform  with  amplitude  modulation  described  by 
this  time  sequence  with  a  direct  linear  interpolation  between  sequence 
points  is  shown  in  Figure  1(c). 

The  waveform  in  figure  1(c)  is  described  by  jr ^  (r)  =  xft )  sin  (10  nr),  so  that 
the  frequency  of  the  modulated  waveform  is  five  times  the  sample  fre¬ 
quency.  The  baseband  time  waveform,  xft),  is  obtained  using  a  linear  inter¬ 
polation  between  sample  points.  The  use  of  waveforms  such  as  sine 
functions  or  Nyquist  pulses  has  been  avoided  because  these  bases  derive 
from  linear  signal  theory  [9].  In  general,  the  best  basis  to  use  to  reconstruct 
a  time-sampled  chaotic  trajectory  is  precisely  the  curve  that  projects  the  tra¬ 
jectory  from  the  current  sample  point  to  the  next  one  on  the  surface  of  sec¬ 
tion.  In  this  case,  the  system  does  not  refer  to  a  continuous  time  system,  so 
the  choice  of  basis  is  arbitrary,  and  I  have  used  a  particularly  simple  one. 
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Fig-jre  1.  (a)  iCO-point  (a) 
sequence  generated  by 
natural  dynamics  of 
map  xMtj  =  Zx„  mod  1, 

(b)  power  spectral 
density,  and 

(c)  modulation 
waveform. 


3.  Forming  a  Pulse 

Because  we  can  describe  any  allowed  dynamical  motion  of  this  system  by 
specifying  the  bits  in  the  binary  fracu  n  for  an  mitral  point  (and  thus  the 
symbolic  future  of  the  system),  I  con.: :::  pulses  of  the  following  form.  As¬ 
sume  that  the  system  is  first  in  the  off  or  approximately  so.  That  is. 
the  initial  point  is  .t0  =  0.000  ...  000b,b...,b-~.  ...  a  long  string  of  zeros  fol¬ 
lowed  by  a  bit  sequence  to  be  specified.  If  the  :r..ual  string  of  zeros  is  long 
enough,  then  the  system  can  be  considered  rr.  be  irxked  to  the  off  state  by 
an  extremely  tiny  control  that  simply  kicks  the  system  back  toward  zero  as 
the  shrift  dynamics  tries  to  move  it  away.  This  would  correspond,  in  a  physi¬ 
cal  system  or  device,  to  locking  the  system  tn  a  static  but  unstable  equilib¬ 
rium  state,  called  an  unstable  fixed  point.  I  also  consider  that  the  system 
eventually  is  returned,  after  emission  of  the  pulse,  to  the  off  state  Thus,  the 
most  general  form  for  a  pulse  of  this  type  is  specified  by  the  initial  point  (or 
symbol  sequence)  x0  =  0.0s  PT 0",  where  0s  denotes  a  precursor  string  of 
zeros  of  length  .V,  0"  represents  an  infinite  tail  string  of  zeros,  and  PT  is  an 
arbitrary  string  of  ones  and  zeros  of  integer  length  T  representing  the  pulse 
bits.  1  assume  that  the  precursor  string  0s  is  arbitrarily  long  for  computa¬ 
tional  purposes,  so  that  the  system  is  very  dose  to  the  off  state  before  the 
pulse  is  initiated. 

The  problem  now  is  simply  to  specify  the  pulse  string  Pr.  The  simplest 
pulse  that  can  exist  occurs  when  PT  =  1;  that  is,  the  pulse  stnng  is  one  bit 
long  and  is  given  by  the  binary  digit  1. 1  will  call  this  the  unitpuhe.  for  ob¬ 
vious  reasons,  but  it  should  not  be  confused  with  the  unit  amplitude  impulse 
in  linear  signal  theory  [9]  This  unit  pulse,  as  will  be  shown,  is  a  growing 
exponential.  It  is  straightforward  to  describe  the  pulse  amplitude  as  a  func¬ 
tion  of  time  (in  this  case  the  integer’  index  r):  Because  the  single  pulse  bit  is 
at  first  deep  in  the  number  r0.  the  system  state  is  initially  off  (practically; 
As  the  long  string  of  zeros  in  the  precursor,  with  the  one  bit  deep  inside,  get 
shifted  out  (to  the  left),  the  system  state  begins  to  move  away  from  off. 
Equivalently,  the  signal  amplitude  builds  from  zero  and  grows  exponen¬ 
tially.  The  shifting  of  this  one  bit  from  deep  inside  the  number  to  succes¬ 
sively  higher  significant  bit  slots  means  that  the  pulse  signal  as  a  function  of 
time  is  given  by  x,  -  \[2  2',  t  £  0;  and  x,  -  0,  t  >  0. 

If  I  define  the  step  function  s,  to  be  unity  for  /  <  0,  and  zero  for  t  >  0,  the  unit 
pulse  can  be  written  as  u,  =  2'~,sl  The  pulse  amplitude  thus  builds  up  expo¬ 
nentially  with  time,  and  then  cuts  off  abruptly  at  t  ~  0.  This  dynamics  is 
portrayed  graphically  in  figure  2,  which  shows  the  graph  of  the  shift  map 

*/  will  use  the  variable  t  to  represent  the  discrete-time  index,  instead  of  a  commonly  used  integer  vanable  like  n.  1  will 
also  use  the  terminology  waveform  when  referring  to  the  sequence  generated  by  the  shift  map.  Both  of  these  uncon¬ 
ventional  usages  were  chosen  to  avoid  referring  to  the  sequence  generated  by  the  shift  as  a  signal,  which  has  the  con¬ 
notation  of  information  transmission  and  signal  processing  This  note  deals  wuh  pulse  generation,  possibly  for  trans¬ 
mission  from  a  high-powered  microwave  source  or  laser,  but  not  intended  for  information  transmission 
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Figure  2.  I  ni t  puhe 
djn.imifs  in  m;i  r  anc. 


along  with  a  map  plane  (function  space)  trajectory  of  the  system.  The  dy¬ 
namics  of  the  state  point  is  represented  by  the  dynamics  cf  a  point  in  the 
map  plane  that  successively  travels  between  the  map  function  and  the  iden¬ 
tity  line  x..,  =  x..  The  reflection  of  the  point  off  the  map  thus  denotes  the  in¬ 
sertion  of  the  value  cf  x,  into  the  map.  Tne  reflection  of  the  point  off  the 
identity  line  represents  the  setting  of  the  value  cf  x,  to  the  value  of  x..-  from 
the  last  iterate.  The  system  point  first  exits  the  off  state  .along  the  dotted 
line,  and  the  amplitude  builds  exponentially  until  it  hits  the  point  .r  = 
1/2  and  locks  b3ck  to  the  off  state.  Note  that  the  system  locks  precisely  to 
off  because  all  the  bits  following  the  pulse  bit  are  zero. 


The  spectrum  of  this  pulse  is  easily  computed: 


LYw)  = 


u.c 


so  that 


v(cj)=  ' 


and  using  the  closed  form  expression  for  the  power  series  yields  t*(  to'  = 
1/(2  -  eu).  The  power  spectrum  is  thus  given  by  }L*(<t»>j2  =  1  (5-4  ecs  <jt. 

The  unit  pulse  waveform  and  its  normalized  power  spectrum  arc  shown  m 
figures  3(a)  and  (b).  respectively.  Note  that  the  power  spectrum  of  this 
time-limited  pulse  is  the  same  as  the  power  spectral  density  for  the  station- 
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Figure  ?  Unit  pulse 

(a)  time  sequence, 

(b)  power  spectrum, 
and  fc)  modulation 
waveform. 
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6  3 


ary  crgodic  signal  produced  by  the  natural  dynamics  of  the  system.  This 
equivalence  provides  a  good  example  that  two  signals  need  not  be  the  same 
to  yield  the  same  power  spectrum,  and  shows  a  deeper  mathematical  con¬ 
nection  between  the  natural  dynamics  and  the  unit  pulse.  (A  natural  trajec¬ 
tory  of  the  system  can  be  described  as  a  random  superposition  of  unit 
pulses;  thus  the  power  spectra  can  be  shown  to  be  the  same.)  An  amplitude 
modulation  waveform  generated  by  this  pulse  is  shown  in  figure  3(c),  given 
by  Uj(r)  ~  u(:)  sin  (10  ,t:),  where  u(r)  is  the  linearly  interpolated  unit  pulse. 

More  complex  pulses  can  be  formed  similarly  by  specifying  the  pulse  string 
PT  Some  pulses  of  length  T  with  different  pulse  strings  are  equivalent.  This 
occurs  because  any  zeros  at  the  beginning  or  end  of  the  pulse  string  do 
nothing  but  shift  the  time  of  the  pulse,  an  inconsequential  detail  in  view  of 
the  fact  that  the  precursor  string  is  arbitrarily  long.  Therefore,  1  adopt  the 
convention  that  all  pulse  strings  must  begin  and  end  in  a  one.  All  bits  be¬ 
tween  the  first  and  last  one  are  then  arbitrary,  and  the  number  of  unique 
pulses  of  length  T  is  given  by  NT  =  2T'2.  Some  of  these  are  symbolic  time- 
reversed  copies  of  another,  but  the  state  point  dynamics  differ  because  of 
the  different  numerical  significance  of  the  bit  positions. 
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4.  Spectral  Shaping 


More  complex  pulses  can  be  described  in  the  context  of  spectral  shaping. 
(There  are  other  ways  to  differentiate  between  the  characteristics  of  the 
pulses,  but  power  spectra  have  become  entrenched  as  a  central  concept  in 
linear  signal  theory. 1  The  power  spectrum  also  has  the  appeal  of  telling 
what  spectral  bandwidth  is  needed  foi  transmission  of  a  given  amount  of 
pulse  energy.  It  is  obvious,  at  least,  that  the  symbolic  description  is  too  de¬ 
tailed  for  differentiating  global  properties  of  pulses.  For  a  long  pulse  length, 
T,  there  are  many  pulses  that  are  nearly  the  same  in  the  sense  that  they  ap¬ 
pear  quasirandom  within  the  pulsewidth. 

During  the  buildup  phase  of  a  pulse,  the  amplitude  grows  exponentially. 
Because  the  pulse  string  is  sequentially  amplified  by  a  factor  of  two,  if  I  let 
the  r  cal  number  A0  =  0 ,PT,  then  the  pulse  amplitude  as  a  function  of  time, 
during  the  buildup  phase,  is  given  by  A(t)~  <0.  At  r  =  0,  the 

leftmost  bit  of  the  pulse  string  is  in  the  most  significant  bit  (msb)  slot  of  the 
binary  fraction  (the  1/2  place).  For  r  >  0,  the  truncation  caused  by  the 
modulo  operation  causes  the  pulse  amplitude  to  change  in  a  chaotic  fashion 
(at  least  for  most  pulse  strings)  until  the  last  bit  (required  to  be  a  1)  is  in  the 
msb  slot.  Then,  the  pulse  amplitude  is  1/2,  and  after  the  next  iteration,  the 
system  again  returns  abruptly  to  the  off  state. 

Now  let  pt  =  •  where  uf  is  a  unit  pulse  that 

reaches  its  maximum  value  of  1/2  at  r  =  k,  and  ak  is  either  zero  or  one  de¬ 
pending  on  the  corresponding  bit  in  the  desired  pulse  string.  Following  this 
procedure,  any  of  the  ZT~ 1  possible  pulses  can  be  constructed  from  the  unit 
pulses,  in  a  manner  roughly  analogous  to  the  formation  of  arbitrary  pulses 
in  linear  signal  theory  from  the  unit  impulse.*  In  this  case,  however,  it  is 
important  to  remember  that  the  superposition  is  occurring  bitwise  in  the  bi¬ 
nary  fraction  symbolic  description  of  the  system  state,  and  that  the  linear 
superposition  always  adds  bits  that  are  never  both  ones.  Thus,  in  the  lan¬ 
guage  of  arithmetic,  no  carries  ever  happen.  If  one  tries  to  use  superposition 
of  two  pulses  that  cause  a  carry,  the  linear  superposition  breaks  down. 

The  possibility  of  using  linear  superposition  like  this  helps  immensely  for 
the  problem  of  spectral  shaping.  Because  all  complex  pulses  are  formed 
from  unit  pulses,  the  spectrum  of  an  arbitrary  pulse  can  be  computed  easily. 


Vi  is  interesting  that  the  unit  impulse  from  classical  signal  theory  is  obtained  in  the  limit  of  a  large  symbol  alphabet 
from  the  concept  introduc'd  here  of  linear  superposition  of  symbol  strings.  If  a  symbol  sequence  0^  d^‘CL^'  ...from 
a  symbol  alphabet  of  cardinality  M  is  represented  in  its  M -ary  fraction  representation,  then  the  msb  slot  containing 
m/  is  M  times  more  significant  than  the  next  most  significant  slot  containing  my.  Thus,  as  M  becomes  eery  large,  >he 
amplitude  of  the  shift,  u,hich  is  given  by  0  m/mimy  .  remains  very  close  to  zero  if  mi  =  0,  and  is  close  to  one  if  m/  = 
M  -  1 .  The  signal  amplitude  is  thus  largely  determined  by  the  msb  slot  of  the  shift,  and  the  unit  pulse  for  the  Si-ary 
alphabet  is  a  sharply  rising  exponential  that  becomes  the  classical  unit  impulse  as  M  — *  •*» . 


Thus 


^-»**«* 


d«* 


Using  the  expression  for  the  unit  pulse,  and  exchanging  the  order  of  the 
sum  yields 


^)=x;rKx;::_2'-^- 


The  sum  over  time  is  just  the  Fourier  transform  of  the  unit  pulse  phase- 
shifted  because  the  peak  does  not  occur  at  l  =  0.  Thus  the  spectrum  becomes 


P(o)  = 


(  \ 


The  expression  for  the  power  spectrum  can  be  computed  by  noting  that  the 
squared  magnitude  of  the  amplitude  spectrum  is  equivalent  to  the  matrix 
expression 


iwr = 


ir4>a 


5-4cos(cd) 


where  <t>  =  [0t/]  =  FFt  is  the  matrix  of  products  of  phase  terms  du  = 
F  =  [/j]  is  a  column  vector  of  phase  terms  ft  =  r44*,  and  a  =  [cj  is  the 
column  vector  of  uni:  pulse  coefficients. 

This  expression  for  the  power  spectrum  yields  an  elegant  solution  to  the 
problem  of  spectral  shaping.  Because  any  of  the  ak  can  be  set  to  zero  or  one 
as  desired  the  spectrum  can  be  varied  over  a  large  range  of  possible  spec¬ 
tra.  Noting  that  the  terms  in  the  numerator  of  |/>(o)]‘  are  the  components  in 
a  Fourier  synthesis  of  a  function,  then  the  problem  reduces  to  the  determi¬ 
nation  of  the  Fourier  coefficients  that  will  yield  a  spectrum  close  to  the  one 
desired.  Recall  that  even  the  coefficients  that  are  variable  can  only  be  zero 
or  one,  so  this  is  not  exactly  like  a  classical  problem  of  undetermined  coef¬ 
ficients. 

Instead  of  solving  the  problem  in  detail  (which  would  obscure  the  basic 
idea),  I  will  give  some  specific  solutions  that  alter  the  spectrum  from  the 
one  occurring  with  the  natural  dynamics.  A  simple  pulse  waveform  con¬ 
structed  from  the  unit  pulse  is  given  by  P,B  =  1111111111,  which  is  a  pulse 
of  10  unit  pulses  in  a  row.  This  pulse  waveform,  along  with  its  power  spec¬ 
trum  and  corresponding  amplitude-modulated  waveform,  is  shown  in  figure 
4.  The  power  density  has  now  become  concentrated  near  o)  =  0.  The  pulse 
therefore  requires  much  less  spectral  bandwidth  for  transmission  than  either 
the  unit  pulse  or  the  signal  produced  by  the  natural  dynamics. 
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Other  pulse  waveforms  can  be  constructed  similarly.  Figure  5  shows  the 
pulse  given  by  P4  =  1010101010,  along  with  its  power  spectrum  and  ampli¬ 
tude  modulation  waveform.  (This  is  a  T  =  9  pu’se  because  the  last  bit  is 
zero.)  Note  that  the  power  spectrum  of  the  pulse  now  has  significant  power 
density  near  cj  -  ±n.  This  is  caused  by  the  presence  of  the  frequency  com¬ 
ponent  near  a  period  two  orbit  of  the  map.  Finally,  figure  6  shows  the  pulse 
P'°  =  1100110011,  along  with  its  power  spectrum  and  amplitude  modula¬ 
tion  waveform.  This  pulse  comes  close  to  a  period  four  orbit,  and  the  spec¬ 
trum  is  enhanced  near  co  =  ±rJ2.  The  bit  patterns  indicated  by  these  pulses, 
if  continued  indefinitely,  would  cause  the  system  to  become  locked  in  a  pe¬ 
riodic  orbit,  period  two  and  four,  respectively.  With  very  long  pulses,  the 
corresponding  Fourier  component  would  become  prominent,  and  would  ap¬ 
proach  a  delta  function  in  the  frequency  domain.  The  laige  lobe  at  co  =  0 
visible  in  all  the  pulse  spectra  in  this  report  would  also  approach  a  delta 
function.  This  lobe,  representing  the  positive  offset  component  of  all  the 
pulses,  would  correspond  to  the  dc  component  of  the  signal  for  long  pulses. 
Thus,  one  can  think  of  the  other  spectral  lines  as  representing  the  amplitude 
modulation  on  the  dc  component  of  the  signal. 
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